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Abstract 

We analyzed mtDNA control region sequences of green turtles (Chelonia mydas) from Arvoredo Island, a foraging 
ground in southern Brazil, and identified eight haplotypes. Of these, CM-A8 (64%) and CM-A5 (22%) were dominant, 
the remainder presenting low frequencies (< 5%). Haplotype (h) and nucleotide (n) diversities were 0.5570 ± 0.0697 
and 0.0021 ± 0.0016, respectively. Exact tests of differentiation and AMOVA <1> ST pairwise values between the study 
area and eight other Atlantic foraging grounds revealed significant differences in most areas, except Ubatuba and 
Rocas/Noronha, in Brazil (p > 0.05). Mixed Stock Analysis, incorporating eleven Atlantic and one Mediterranean 
rookery as possible sources of individuals, indicated Ascension and Aves islands as the main contributing stocks to 
the Arvoredo aggregation (68.01% and 22.96%, respectively). These results demonstrate the extensive relation- 
ships between Arvoredo Island and other Atlantic foraging and breeding areas. Such an understanding provides a 
framework for establishing adequate management and conservation strategies for this endangered species. 
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The green turtle (Chelonia mydas) is a marine reptile 
of worldwide tropical and subtropical distribution, cur- 
rently classified by the World Conservation Unit as endan- 
gered (IUCN, 2007). These animals present complex and 
long life histories, together with a highly migratory behav- 
iour (Meylan, 1995; Godley et al, 2003). Due to the large 
temporal and spatial scales involved, various aspects of 
their life cycle are quite difficult to elucidate by conven- 
tional approaches, and must be solved by using indirect re- 
search methods, such as molecular genetics (Avise, 2007; 
Bowen and Karl, 2007). 

Mitochondrial DNA (mtDNA) control region studies 
have been increasingly applied to marine turtles, whereby 
the development of genetic tags for these animals has con- 
tributed to the acquisition of valuable data on their molecu- 
lar evolution, population structure, reproductive behavior 
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and migration ecology, besides providing a foundation for 
conservation and management strategies (Moritz, 1994; 
Avise, 2007; Bowen and Karl, 2007). In this context, green 
turtles have emerged as model organisms for such studies 
(Avise, 2007). These animals forage in "mixed stocks" 
composed of individuals from several cohorts and from 
various nesting beaches (rookeries) which aggregate at 
feeding grounds (Bass and Witzell, 2000; Bass et al, 2006; 
Avise, 2007; Bowen and Karl, 2007). Due to the phylo- 
patric behaviour of the females of this species, nesting 
assemblages are genetically structured in terms of mater- 
nally-inherited characters, thereby permitting the evalua- 
tion of the natal origins of individuals found in mixed ag- 
gregations (Bowen, 1995; Bowen and Karl, 2007). The 
assessment of the genetic composition of mixed aggrega- 
tions is currently one of the research priorities for this spe- 
cies (Formia et al, 2006). This, together with the determi- 
nation of the relationships among foraging and breeding 
populations of sea turtles, are essential for constituting se- 
cure guide lines in the development of successful conserva- 
tion strategies for these endangered animals. 
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Sampling was undertaken at Arvoredo Island, located 
within the Arvoredo Marine Biological Reserve (27° 17' S 
and 48° 22' W), in July 2005, January-February 2006 and 
July 2006, at five different sites located on the western and 
northern parts of the island (Figure 1). Tissue samples were 
obtained from the flippers of 49 juvenile green turtles 
hand-captured through free and SCUBA dives, by using 
5 mm disposable biopsy punches. The samples were then 
conserved in absolute ethanol and kept at -20 °C. Curved 
carapace length and the weight of sampled individuals 
ranged from 35 to 72.5 cm (mean 52 cm) and 7.5 to 45 kg 
(mean 19.9 kg), respectively. 

DNA extraction was performed through the standard 
phenol xhlorophorm method with precipitation in absolute 
ethanol (Hillis etal, 1996). Control region fragments were 
amplified via polymerase chain reactions (PCR) using the 
primers LTCM1 and HDCM1 (Allard et al, 1994), under 
the following conditions: initial denaturation of 1 min at 
94 °C, 35 cycles of 30 s at 94 °C, 1 min at 50 °C, 1 min at 
72 °C, and a final 5 min extension at 72 °C. Products were 
purified using Illustra GFX purification kits (GE Health- 
care, U.S.A.), and sequenced in both directions using 
DYEnamic ET dye terminator kits in a MegaBACE 
500 DNA sequencer (GE Healthcare, U.S.A.). 

Sequences (491 bp) were aligned by means of Clustal 
X 1.83 software (Thompson et al, 1997), and haplotypes 
classified according to the Archie Carr Center for Sea Tur- 




Figure 1 - Location of Arvoredo Island (Al - triangle) and other foraging 
and nesting areas used for comparison and Mixed Stock Analysis. Abbre- 
viations for foraging grounds (squares) are: UB (Ubatuba), R/N (Ro- 
cas/Noronha), AF (Almofala), BA (Barbados), BH (Bahamas), NI (Nica- 
ragua), FL (Florida) and NC (North Carolina). Abbreviations for nesting 
areas (circles) are: TI (Trindade Island), R/N (Rocas/Noronha), AS (As- 
cension Island), GB (Guinea Bissau), BI (Bioko), ST (Sao Tome), AV 
(Aves Island), SU (Surinam), MX (Mexico), CR (Costa Rica), FL (Flo- 
rida) and CY (Cyprus). 



tie Research online genetic bank (Florida University). A 
minimum spanning network demonstrating relationships 
among haplotypes was set up using TCS 1.3 software 
(Clement et al, 2000). Exact tests of differentiation be- 
tween Arvoredo Island and other Atlantic foraging grounds 
were carried out with Arlequin 3.11 (Excoffier et al. , 2005), 
using Markov Chain Monte Carlo (MCMC) of 10000 per- 
mutations with 1000 dememorization steps. Pairwise <E>- 
statistics (<D S t, which summarizes the degree of differentia- 
tion between populations) were checked through Analysis 
of Molecular Variance (AMOVA) conducted with 10000 
permutations with Arlequin 3.11, according to the Tamu- 
ra-Nei model of nucleotide substitution. The Brazilian for- 
aging grounds included in these analyses for comparison 
were Ubatuba (SP), Almofala (CE) (Naro-Maciel et al, 
2007), Rocas Atoll (RN) and Fernando de Noronha (PE) 
(Bjorndal et al, 2006). The latter two were grouped into 
one single unit for all analyses, due to geographic proximity 
(c.a. 150 km) and small sample size, being hereafter re- 
ferred to as Rocas/Noronha. Nicaragua (Bass et al, 1998), 
Barbados (Luke et al, 2004), Bahamas (Lahanas et al, 
1998), Florida (Bass and Witzell 2000) and North Carolina 
(Bass et al, 2006), in the Caribbean and North Atlantic, 
were also included for comparison. Structuring between 
foraging areas grouped into North and South Atlantic ag- 
gregations was checked through AMOVA. 

Probable natal origins were defined through Mixed 
Stock Analysis (MSA) employing Bayes software (Pella 
and Masuda, 2001), and considering equal prior probabili- 
ties assigned to each source. Source populations employed 
as possible contributors to the Arvoredo Island group corre- 
spond to all the Atlantic and Mediterranean rookeries de- 
scribed in literature by Bjorndal et al. (2005, 2006), Formia 
et al. (2006, 2007), Encalada et al. (1996) and Kaska 
(2000), namely, Trindade Island and Rocas/Noronha (Bra- 
zil), Ascension Island (United Kingdom), Poilao (Guinea 
Bissau), Bioko Island (Equatorial Guinea), Sao Tome (De- 
mocratic Republic of Sao Tome and Principe), Aves Island 
(Venezuela), Matapica (Surinam), Quintana Roo (Mexico), 
Tortuguero (Costa Rica), Florida (United States) and Lara 
Bay (Cyprus). Principe (Democratic Republic of Sao Tome 
and Principe) was excluded from this analysis due to the 
small size of the sample. 

We encountered eight polymorphic sites which de- 
fined eight previously described Atlantic Ocean haplo- 
types. The predominant haplotype was CM-A8 (64%), 
commonly found in South Atlantic rookeries, followed by 
CM-A5 (22%), mainly found in the Costa Rica, Surinam 
and Aves Island rookeries. The remaining haplotypes were 
relatively rare, with less than 5% frequency. To date, CM- 
A9 (2%), CM-A24 (4%) and CM-A32 (2%) have only been 
observed in the South Atlantic rookeries of Rocas Atoll, 
Trindade and Ascension Island, whereas CM-A10 (2%) has 
been found in Rocas Atoll and Ascension Island. CM-A39 
(2%), previously unregistered in foraging areas, and CM- 
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ing areas, with the exception of Ubatuba and Rocas/No- 
ronha in Brazil (p = 0.4776 and 0.3077, respectively). Simi- 
lar results were inferred from AMOVA (p = 0.1292 and 
0.6261). By grouping foraging aggregations into North and 
South Atlantic and using AMOVA, strong structuring was 
revealed (<J> ST = 0.6913 p < 0.01). From MSA, it was indi- 
cated that Ascension and Aves Islands are the main contri- 
butors to the Arvoredo aggregation, with lesser contribu- 
tions from Guinea Bissau and Trindade Island, as shown in 
Table 2. 

High CM-A8 frequency in the study area is in accor- 
dance with the predominance of this haplotype in various 
nesting and feeding areas in the Atlantic, and is consistent 
with the suggestion of it being the closest relative to an an- 
cestral haplotype in the Atlantic basin. Haplotype CM-A5 
was the second most frequent, as was noted in other south 
Atlantic feeding grounds, and in accordance with its high 
frequency in large Caribbean rookeries (Bjorndal et ah, 
2005, 2006; Formia et ah, 2006, 2007; Naro-Maciel et ah, 
2007). Elevated h values are found in most green turtle for- 
aging areas, as expected when considering that these aggre- 
gations are composed of mixed stocks (Bass and Witzell, 
2000). Low 7i values were also expected due to the slight 
variation observed between haplotypes. 

The distribution of haplotypes among foraging 
grounds is apparently non-random, with significant differ- 
entiation among individual areas and strong structuring 
between North and South Atlantic aggregations. The life 
history patterns of sea turtles may account for such structur- 
ing, with the pelagic stage determining the areas to which 
these animals will recruit, possibly at the whim of ocean 
currents (Musick and Limpus, 1997; Luschi et ah, 2003). 
Arvoredo Island was not significantly different from the 
closest genetically-described southwestern Atlantic forag- 
ing ground, Ubatuba (ca. 755 km), thereby indicating that 
foraging areas can present similarity in mtDNA at small 
spatial scales. Such a hypothesis is corroborated by Al- 

Table 1 - Haplotype (h) and nucleotide (n) diversity estimates + standard deviations for all compared foraging aggregations. 



Foraging ground 


Haplotypes 


h 




Sample size 


Arvoredo Island 


8 


0.5570 + 0.0697 


0.0021 ±0.0016 


49 


Ubatuba" 


10 


0.4460 ± 0.0556 


0.0020 ±0.0015 


113 


Rocas/Noronha b 


6 


0.5887 + 0.0911 


0.0019 ±0.0015 


32 


Almofala" 


13 


0.7168 ±0.0306 


0.0067 ± 0.0039 


1 17 


Barbados 0 


8 


0.7734 ±0.0276 


0.0105 ±0.0057 


60 


Bahamas' 1 


6 


0.3703 ±0.0650 


0.0066 ± 0.0038 


79 


Nicaragua 6 


2 


0.1831 ±0.0621 


0.0039 ± 0.0025 


60 


Florida f 


6 


0.4855 ± 0.0668 


0.0032 ±0.0021 


62 


North Carolina 8 


8 


0.6778 ±0.0310 


0.0052 ±0.0031 


106 


Average 


7 


0.5334 


0.0040 


68 



"Naro-Maciel et al. 2007. b Bjorndal et al. 2006. c Luke et al. 1994. d Lahanas et at. 1998. e Bass et al. 1998. f Bass and Witzell 2000. E Bass et al. 2006. 



A45 (2%), with only one register m feeding grounds, have 
been described only in animals from the Ascension Island 
rookery. Haplotypes were distinguished by a maximum of 
two variations, as shown in the Minimum Spanning Net- 
work (Figure 2). 

Haplotype (h) and nucleotide (jx) diversity estimates 
encountered for the study area were 0.5570 + 0.0697 and 
0.0021 + 0.0016, respectively. Diversity estimates for Ar- 
voredo Island and other Atlantic foraging grounds are listed 
in Table 1 . Exact tests of differentiation based on haplotype 
frequencies demonstrated general differentiation among all 
feeding areas (p = 0.000). According to these tests, Ar- 
voredo Island was significantly different from most forag- 




Figure 2 - Minimum spanning network of mtDNA control region relation- 
ships encountered at Arvoredo Island. Hash lines represent 1 basepair sub- 
stitution between haplotypes. 
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Table 2 - Mixed stock analysis based on Bayesian methods considering equal priors, with mean, standard deviation (S.D.), 2.5% quantile, median and 
97.5% quantile. 



Stock 


Mean 


S.D. 


2.5% 


Median 


97.5% 


Trindade Island" 


0.0218 


0.0535 


0.0000 


0.0001 


0.1852 


Rocas/Noronha" 


0.0161 


0.0471 


0.0000 


0.0000 


0.1700 


Ascension Island 11, c- d 


0.6801 


0.1171 


0.3869 


0.7029 


0.8407 


Guinea Bissau" 


0.0197 


0.0542 


0.0000 


0.0000 


0.1948 


Bioko c 


0.0174 


0.0504 


0.0000 


0.0000 


0.1710 


Sao Tome 0 


0.0062 


0.0220 


0.0000 


0.0000 


0.0663 


Aves Island* 1 


0.2296 


0.0597 


0.1257 


0.2257 


0.3592 


Surinam* 1 


0.0019 


0.0064 


0.0000 


0.0000 


0.0199 


Mexico d 


0.0019 


0.0064 


0.0000 


0.0000 


0.0196 


Costa Rica" 


0.0019 


0.0063 


0.0000 


0.0000 


0.0193 


Florida* 1 


0.0017 


0.0058 


0.0000 


0.0000 


0.0177 


Cyprus df 


0.0017 


0.0056 


0.0000 


0.0000 


0.0159 



"Bjorndal et al. 2006. b Formia et al. 2007. Tormia et al. 2006. d Encalada et al. 1996. e Bjorndal et al. 2005. f Kaska 2000. 



mofala, the most distant southwestern Atlantic foraging 
ground from Arvoredo Island (ca. 3800 km), being signifi- 
cantly different from the study area. This difference could 
also be due to its proximity to the Caribbean region, with its 
strong structuring within the Atlantic Ocean (Bass et al. 
2006). The similarity observed between relatively close 
feeding grounds could possibly be attributed to movements 
along the coast, which may be influenced by factors such as 
current intensity, variations in temperature, food availabil- 
ity and continuous recruitment to coastal zones (Bass et al, 
2006). Despite many animals presenting high fidelity to 
foraging areas, it is known that non-reproductive costal 
movements of juvenile green turtles may occur (Godley et 
al, 2003; Bass et al, 2006), the geographic nearness of the 
areas and major coastal currents also possibly constituting 
important factors in these movements. Green turtles present 
at least short-term fidelity to Arvoredo Island, as demon- 
strated by various recaptures over a three-year study period 
(Reisser et al, 2008). Nevertheless, one animal tagged in 
the area was encountered six months later by members of 
Project Tamar-ICMBio, stranded at Caraguatatuba in Sao 
Paulo state, over 700 km away, thus giving evidence of 
non-reproductive migration in coastal waters. Migration 
between Sao Paulo and southern Brazil has also been ob- 
served by Marcovaldi et al (2000), in which a green turtle, 
initially tagged at Ubatuba, was recaptured three months 
later in Bombinhas, SC, only 1 0 km from Arvoredo Island. 

As was the case for other south Atlantic foraging ar- 
eas (Bjorndal et al, 2006; Naro- Maciel et al, 2007), the 
main stock contributing to the Arvoredo aggregation be- 
longs to Ascension Island. Green turtle movements be- 
tween Ascension and Brazil have often been noted through 
mark-recapture and telemetry studies (Meylan, 1995, Lus- 
chi et al, 1998, Hays et al, 2002). The large nesting popu- 
lation (the second largest in the Atlantic, with 



approximately 3800 females nesting annually; Broderick et 
al. , 2006) and favorable ocean currents are the most proba- 
ble explanations for such a high contribution. The second 
largest contributor was Aves Island, although there is a lack 
of tagging evidence on migrations to-and-from Brazil. 
However, the relatively large rookery size (300-500 fe- 
males nesting annually; Seminoff, 2002), and the strong 
link between Caribbean rookeries and Brazilian foraging 
grounds, as demonstrated by Lima et al. , (2008), give sup- 
port to this conclusion. The connection between African 
rookeries and Brazilian foraging grounds is not evident, 
possibly due to the limited number of studies dealing with 
the African continent. Estimates inferred from MSA indi- 
cated that African contributions as a whole to Arvoredo Is- 
land were generally low, although those from Guinea 
Bissau and Bioko were relatively high compared to other 
African nesting areas. Naro-Maciel et al. (2007) also ob- 
served a relatively high contribution from Guinea Bissau to 
Ubatuba. This could be a consequence of the fixed charac- 
teristics of this area for the commonly found haplotype 
CM-08 (Formia et al, 2006), which could have affected 
MSA estimates. Bioko also presents a high frequency of 
haplotype CM-08 (90%), also possibly interfering with the 
analysis. The contribution from Trindade Island is appar- 
ently underestimated when considering that this island sup- 
ports the largest nesting area in Brazil (approximately 
300-400 females during the last nesting season - Soares LS, 
personal communication to PLR), and is the nearest rook- 
ery to the study area (ca. 2100 km). Furthermore, numerous 
recaptures of green turtles tagged in this area have been reg- 
istered along the Brazilian coast (Marcovaldi et al, 2000). 
Low estimated contributions from Trindade Island have 
also been registered for the previously cited mixed stocks 
described in Brazil (Almofala, Ubatuba, Rocas/Noronha). 
However, in a recent study by Bolker et al. (2007), a 
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'many-to-many' MSA approach with the incorporation of 
multiple mixed stocks gave evidence of higher contribu- 
tions from Trindade Island to northeastern Brazil than those 
previously published. This could corroborate the hypothe- 
sis that Trindade 's contribution to the study area is underes- 
timated. Nonetheless, further investigation is necessary to 
clarify this. 

The assumption that all sources and all mixtures are 
well described is a great problem with MSA, since this is not 
always the case. The presence of foraging ground haplotypes 
which have not been described at nesting areas clearly indi- 
cates that some rookeries may be inadequately described or 
not even at all, as was noted by Bass et al. (2006), Formia et 
al. (2007) and Naro-Maciel et al. (2007). Furthermore, ha- 
plotypes being encountered in rookeries but not in foraging 
areas demonstrates insufficient research at feeding grounds. 
Therefore, this analysis should be interpreted together with 
all available evidence {i.e. demographic, ecological, and mo- 
lecular), in order to reach conclusive information on the life 
history patterns of sea turtles. 

Describing the genetic characteristics of juvenile 
green turtle foraging grounds and defining their relation- 
ship with other feeding and breeding grounds provide a 
framework for successfully conserving and managing this 
species. The extensive Brazilian coastline and oceanic is- 
lands harbor countless foraging grounds, besides three 
rookeries of which two are relatively large, thereby urging 
investigation and protection for conservation purposes. Im- 
pacts affecting foraging areas may also influence distant 
rookeries. Thus, the protection of feeding zones could be a 
big step towards the protection of their contributing stocks. 
The distribution and migrations of green turtles surpass na- 
tional boundaries, wherefore urging national and interna- 
tional efforts and cooperation is essential for assuring the 
survival of this species. 
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